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Abstract 

We measure spin transport in high mobility suspended graphene {pL « 10 5 cm 2 /Vs), obtain- 
ing a (spin) diffusion coefficient of 0.1 m 2 /s and giving a lower bound on the spin relaxation 
time (t v « 150 ps) and spin relaxation length (A s =4.7 /im) for intrinsic graphene. We develop 
a theoretical model considering the different graphene regions of our devices that explains our 
experimental data. 

The prospectives for graphene spintronics are very positive, with theoretical predictions of 
spin relaxation times (t s ) in the order of hundreds of nanoseconds and higher 1 2 due to the lack of 
nuclear spin for carbon's most common isotope ( 12 C) and weak intrinsic spin-orbit coupling. On 
the other hand the experimental results typically find x s of hundreds of picoseconds J^l In order to 
clarify the limitations and mechanisms for spin relaxation in monolayer graphene devices a lot of 
effort has been done by both experimentalists^! a nd theoreticians' 2 ^^ but up to now this topic 
is still under debate. Theoretical results pointed out that the D'Yakonov-Perel mechanism^ and 
*To whom correspondence should be addressed 
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locally enhanced spin-orbit coupling due to impurities and lattice deformations' 2 ^* are the most 
probable mechanisms for spin relaxation for measurements in graphene spintronic devices.^ * 11 ' 12 ' 
Since all devices studied so far were fabricated on substrates (S iO#! or SidH) in which the 
intrinsic properties of graphene are masked due to the electronically coupling to the substrate,^ 
no previous study was able to confirm these recent theoretical predictions. In order to address 
this problem we study spin transport in high-mobility suspended graphene devices. By removing 
the substrate, thereby suspending the graphene flake, we are not only capable of achieving a high 
quality device with low contamination,^^ but we can also investigate how the absence of the 
rough SiC>2 substrate influences the spin transport in graphene. Moreover it opens the possibility 
to exploit the exquisite mechanical properties of graphene^ and to study how pseudo-magnetic 
fields and strainL^2l a ff ec t spins in graphene, paving the way for graphene in the field of spin- 
nanomechanical applications.^^ 

The effective spin-orbit (SO) field that the spins experience in graphene is directly related to 
the density of impurities and adatoms.^^This effective SO field causes a momentum-dependent 
spin precession, resulting in spin decoherence and relaxation. This mechanism is known as the 
D'Yakonov-Perel mechanism for spin relaxation. 13 Assuming such a mechanism for spin relax- 
ation 24 with an effective SO coupling A$o, the spin relaxation time behaves like:* 12 * 13 * 

1 4A 2 

- = -#T P (1) 

where x p is the momentum relaxation time and h the reduced Planck constant. With the reduction 
of adatoms and impurities we should obtain a low value of A$o approaching the theoretical pre- 
dictions^* 25 * 26 ' (A f ^ ory « 10 jUeV) and a high value for the momentum relaxation time x p > 0. 1 ps. 
Since v s increases quadratically with the reduction in Aso, we expect that for clean samples we 
reach the initial theoretical limits of x s > 10 ns. 

At a first glance, suspending the graphene flake to obtain a very high quality device seems the 
best approach. But the main challenge is that the most common technique to suspend graphene 
flakes™ 1 is acid-based and therefore not compatible with ferromagnetic metals necessary for 



spin transport. The reason is that the acid used to etch the SiC>2 underneath the graphene also 
etches away the ferromagnetic metals used for the electrodes. Also, devices produced by the 
standard technique are typically short, in the order of 1 /im or less to avoid that the graphene flakes 
collapses when a gate-voltage is applied. This is undesirable for spin precession measurements 
since the time for one period of precession for low magnetic fields (B < 1 T) have to be in the same 
order as the diffusion time of the spins in order to obtain a Hanle precession curve^ showing a 
complete spin precession. To overcome these issues we developed a polymer-based method^ in 
which we are able to produce flakes suspended over long distances and contacted by ferromagnetic 
electrodes with highly resistive barriers. In our process the graphene flakes are exfoliated on top 
of a 1 /im thick Lift-Off Resist (LOR) film spin-coated on a Si/Si02 (500 nm) substrate. Single 
layer graphene flakes to be used in our devices are selected by optical contrast.^21221 j^e highly 
doped Si substrate is used to apply a back-gate voltage V g that induces a charge carrier density 
according to n = Oc(V g — Vq), where a = 0.45 x 10 10 cm _2 V _1 is the effective gate capacitance 
and Vb is the position in gate voltage of the minimum of conductance. Using standard electron 
beam lithography (EBL) and metal evaporation methods we deposit the 1 nm thick AI2O3 tunnel 
barriers and the 60 nm Co contacts (Figure la). A second EBL step is used to make parts of 
the graphene flake suspended (Figure 1 b). The detailed device fabrication is described in the 
supporting information. The resistance of our contacts (usually R^ ~20 k£l) and the gate-voltage 
dependence of our devices are characterized at the beginning of every set of measurements. Before 
the quality improvement via the current annealing procedure (described in the next paragraph), the 
devices typically show high p-doping and only a small change of the resistance as a function of 
is observedJ 28 * 31 ! 

To improve the quality of our suspended graphene devices we perform current annealing^ 
in vacuum (at a pressure better than 10~ 6 mbar) at a base temperature of 4.2 K (see supporting 
information). Since the contact resistance of the spin injector and detector in our samples has 
to be kept at high values (>10 to avoid the impedance mismatch and contact induced spin 
relaxation,! 4 * 7 ! 33 ! we use the outer electrodes to apply the high current densities capable of heating 
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the graphene flake to high temperatures (Figure lc). As shown in Figure lb and c, our devices are 
composed of two supported outer parts with a central suspended part. During the annealing the 
graphene flake is cooled by the contacts and the substrate in the supported parts and mainly heats 
up in the suspended region. When a sufficiently high temperature is achieved, the polymer residues 
and contaminants are removed from the suspended central region. Finally the high quality of the 
devices can be verified via the behaviour of the square resistance of the central region (RJ as a 
function of the gate- voltage (Figure le). It is worth noting that the value and the gate dependence of 
the sheet resistance of the supported regions do not change significantly after the current annealing 
procedure. This means that the values of carrier density presented here are representative for the 
central suspended region. We obtained two high quality devices in different flakes with mobilities 
higher than /i = 10 5 cm 2 /Vs. One of the devices was current annealed twice and showed a mobility 
of \l ~ 1.2 x 10 5 cm 2 /Vs after the first current annealing and /i ps 3 x 10 5 cm 2 /Vs after the 
second. This device was also characterized for charge and spin transport in between the two current 
annealing steps. Here we show the results for this representative high mobility sample after the 
second current annealing. The separation between the inner electrodes (suspended region) is L = 
2.5 jLim and the electronic mobility is /i w 3 x 10 5 cm 2 /Vs (Figure le). Similar results for both 
charge and spin transport were obtained in the other devices. 

The values for mobility were obtained in two different ways: either by using the equation 
pL = — and taking the value of /i at a carrier density of n = 1 x 10 10 cm~ 2 , or by fitting the 
conductivity curves with the formula 1/c = l/(ne/l + 0b) + pj 34 l 35 l where o is the graphene 
conductivity, e the elementary charge, Oq the conductivity at the charge neutrality point and p s 
is the contribution of short range scattering. The values obtained by both ways are consistent 
with each other and are in the order of 10 5 cm 2 /Vs. We often observe a background resistance 
in our R sq (V g ) curves (Figure le). The most probable source for this background resistance is the 
non-cleaned region underneath and/or close by the contacts. We did not subtract this background 
resistance in our calculations since this would only lead to higher mobility values. 

To perform the spin-transport measurements we use a non-local geometry where the charge 
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Figure 1 : (a) The graphene flakes are exfoliated on a LOR film spin-coated on a Si/Si02 substrate 
followed by the deposition of Co contacts with aluminium oxide insulating barriers to avoid the 
conductivity mismatch, (b) Local 4-probe measurement geometry to characterize the graphene 
resistance as a function of the gate-voltage applied to the highly doped Si substrate. A second 
electron beam lithography step is performed to suspend the central region, (c) Scheme of the 
current annealing setup used to clean the graphene devices. The high DC bias is applied to the 
outer electrodes to avoid the degradation of the inner contacts since these are to be used for spin 
injection/detection, (d) Scanning electron micrographs of a typical device. The scale bars are 
1 jUm. (e) Local 4-probe resistance (top) and mobility (bottom) versus carrier density and gate- 
voltage after the current annealing procedure at 4.2 K. The carrier density shown corresponds to 
the carrier density in the central suspended region (see main text). 
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current path is separated from the voltage contacts to exclude spurious signals 6 (Figure 2a). When 
we sweep the in-plane magnetic field we can align the two inner ferromagnetic contacts in a par- 
allel or anti-parallel configuration. To extract the spin diffusion coefficient D s and spin relaxation 
time T s , we perform Hanle precession measurements and use the solutions for the Bloch equations 
in the diffusive regime.^^ The Hanle precession measurements are done by measuring the non- 
local resistance as a function of an applied perpendicular magnetic field. To eliminate background 
signals in our analysis we fit the data for the total spin-signal given by: R s = "' 2 , where R n[ ( ' 
is the non-local resistance in the parallel (anti-parallel) configuration of the inner electrodes (Fig- 
ure 2b). The values for D s and v s obtained for the suspended non-annealed samples are typically 
D s « 0.02 m 2 /s and z s « 200 ps, giving a spin relaxation length of X s = \fD s x s ~ 2 /im, showing 
no dependence on the temperature (from room temperature down to 4.2 K). Since the results for 
the non-annealed samples show no substantial difference from measurements done in fully LOR 
supported devices (see supporting information), we can conclude that the substrate is not the main 
factor limiting the spin relaxation for our samples before the cleaning procedure. When we con- 
sider roughness effects, this invariance of the spin relaxation time in our measurements without the 
rough substrate is in agreement with the results by Avsar et al., 5 where it is shown that ripples in 
the graphene flakes have minor (or no) effects on the spin transport parameters. Also, the widths 
of our contacts in our samples are much smaller than the spin relaxation length. So we believe that 
the regions underneath the contacts do not change significantly the spin transport properties in our 
measurements. 

After cleaning by current annealing, we extract T s and D s as a function of the carrier density via 
Hanle precession measurements (Figure 2 c). Comparing our results for the clean high-mobility 
samples to previous studies on SiC>2 supported devices'^ we observe an approximately 3 times 
higher spin diffusion coefficient while the values for the spin relaxation time remain similar. Note 
that the experimental conditions allow us to extract a lower bound for the spin relaxation time in 
the suspended graphene region, this leads to important conclusions as we will discuss below. All 
the spin transport measurements presented here were performed at low temperatures (4.2 K). 
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Figure 2: (a) The 4-probe non-local device configuration for Hanle precession and spin-valve 
measurements. For the precession measurements the two inner contacts are aligned in a parallel 
(anti-parallel) configuration and the non-local resistance is measured as a function of a perpendic- 
ular magnetic field B^. (b) A measured Hanle precession curve (black circles) with a fit using 
the solutions for the Bloch equations (red line). The spin signal R s is recorded as a function of 
the perpendicular magnetic field. Left inset: the non-local resistance as a function of the parallel 
magnetic field. The positive (negative) values of R n i are due to parallel (anti-parallel) alignment 
of the inner electrodes. Right inset: the non-local resistance as a function of the perpendicular 
magnetic field for the parallel (black) and anti -parallel (red) configuration of the contacts, (c) Spin 
and charge diffusion coefficients D s (red squares) and D c (dashed line) respectively, spin relaxation 
time T s and spin relaxation length k s as a function of the carrier density in the central suspended 
region. All the measurements shown were performed at 4.2 K. 
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We start by discussing the results for the spin diffusion coefficient. In Figure 2c we have a 
comparison of the values of D s extracted from the Hanle precession measurements (red squares) 
with the value of the charge diffusion coefficient (dashed line) given by the Einstein relation: 
D c = 1/ (R sq e 2 v(E)), where v(E) is the density of states of graphene at the energy E. In this case 
the square resistance, R sq , was extracted by local 4-probe measurements, as depicted in Figure lb. 
It can be seen that D 4 and D c are in reasonable agreement, which is in accordance to previous 
works on exfoliated monolayer graphene. 3 

We now turn our attention to the results for the spin relaxation time X s (Figure 2c). If the spin 
relaxation time in graphene is limited by a locally enhanced SO coupling^^ due to impurities 
and adatoms, we would expect a high x s for a high mobility sample, since an increased mobility 
is related to a reduction in the density of impurities. Intriguingly, in our results for high mobility 
graphene spin-valves, in which the mean free path is in the order of a micrometer, x s is still in the 
order of hundreds of ps. When we calculate the spin relaxation length A = \fD/z s we obtain large 
values, up to A = 4.7 /im, even with the low spin relaxation time in our samples. But one question 
remains: what limits the measured spin relaxation time in our high-quality devices? To address 
this issue we performed numerical simulations of spin transport in our devices. 

To properly represent our devices we extend the model adopted by Popinciuc et al. 7 We con- 
sider our devices as a three-part system separated by two boundaries. The central part represents 
the suspended region and the two identical outer parts represent the supported regions (Figure 3a). 
The spins are injected at the left boundary by an injector of polarization P c \ and contact resis- 
tance R c \ and detected by a contact at the right boundary with polarization and contact resistance 
P c 2 and R c 2 respectively. The spin diffusion coefficients, spin relaxation times and the conduc- 
tivities, denoted by D^ y T ;(o ) and o { f \ respectively, for the inner (outer) parts can be defined 
independently. We assume that the spin accumulation p. s (x) obeys the Bloch equation for diffu- 
sion in one-dimension with an applied magnetic field B: D s V 2 fl s — & + yB x pi s = 0, where y is 
the gyromagnetic ratio with y = g/ig/i -1 where g the g-factor, /ig the Bohr magneton and h the 
reduced Planck constant. For the boundary conditions we assume that p. s is continuous and goes 
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to zero at x = ±°°, the spin current at the right boundary is continuous and at the left boundary 
it is discontinuous by a value determined by the spin injection. We also include the contact in- 
duced spin relaxation due to back-diffusion,^ although for contact resistances in the order of those 
we encounter in our experiments, our simulations showed no substantial difference from the re- 
sults considering the limit of infinite contact resistances. Details of the model and simulation are 
included in the supporting information. 

By adding a perpendicular magnetic field B = Bz the spins can precess and the voltage at the 
detector electrode is calculated as a function of B, simulating the Hanle precession measurements. 
The obtained data is then fitted the same way we fit our experimental data, with the solutions to 
the Bloch equation for a homogeneous system, and effective values for D s and T s are obtained and 
now denoted as D f it and Xfu respectively. These effective values are then compared to the values 
used in the simulation. 

Since we have four different parameters to consider (T;, t , A and D a ), we change them one by 
one, keeping the others at a constant value. We keep the sheet resistance of the three regions at the 
same value of 1 kQ. The effect of varying the sheet resistance of the inner and outer regions are 
presented in the supporting information. Keeping T, = % = 200 ps and D, = 0. 1 m /s, we varied 
D a from 0.01 to 0.5 m 2 /s and extract Dtu (Figure 3b, blue circles) which we compare to D, (dotted 
black line) and D a (dashed grey line). It can be seen that a change in the spin-diffusion coefficient 
of the outer parts does not influence the results obtained by the Hanle precession fits in the studied 
range. For all values of D () we always obtained a value close to the value of D,, Dfu ^0.1 m 2 /s. 
Keeping now D a constant and varying Dj (Figure 3b, red stars) we confirm that from the Hanle 
precession analysis we always obtain a value Dfj t very close to the actual value of D ; (dashed 
grey line). Such a result agrees with our experimental measurements, where the values obtained 
for D s from the Hanle precession measurements are in good agreement with the charge diffusion 
coefficient, D c , for the inner suspended part. In other words, the diffusion coefficient we measure 
for our inhomogeneous devices is determined by the one in the central high mobility region. 

Having confirmed that the results obtained for the diffusion constant from our experiments and 
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Figure 3: (a) A cartoon of the system used in our simulations: two identical and semi-infinite parts 
connected by one inner region with length L, all of width W. The spins are injected at the left 
boundary (x=0) and detected at the right boundary (x=L). (b) A graph showing the results obtained 
for the diffusion constant by fitting the simulated data with D ; (D ) constant at the value 
represented by the dotted black line and changing D (D,) following the dashed grey line in blue 
circles (red stars), (c) Results for the spin relaxation time Tp t by fitting the simulation data with T ; 
(t ) constant at the value represented by the dotted black line and changing x (t,-) following the 
dashed gray line in blue circles (red stars). The resistivities of the tree regions were kept at 1 \SL 
for the results in (b) and (c). 
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simulations follow the same trend, we now analyze the effect of the spin relaxation times. We fix 
the spin diffusion coefficients at D, = 0.1 m 2 /s and D () = 0.01 m 2 /s, which are approximately the 
values we find in our experiments, and keep T ; - = 200 ps constant (dotted black line). We then vary 
x from 20 to 2000 ps (Figure 3c, blue circles). The values obtained from the fit of the simulated 
data for Xfu follow the increase in x (dashed grey line) and only start deviating at large values 
x ps 2000 ps. On the other hand, when x is kept constant at 200 ps and we change T,, the values 
obtained from the fits seem to be determined by the value of the spin relaxation time in the outer 
regions. From our calculations, this result holds for D„ < Du which falls in our experimental range. 

The results for the spin relaxation time in our model can be explained in a qualitative manner. 
Due to its diffusive motion, a part of the injected spins spend some time in the supported region 
before crossing the suspended region reaching the detector electrode (right boundary). For a suffi- 
ciently low diffusion coefficient in the supported part, the time spent there is enough to relax most 
of the spins in the case of a short spin lifetime. Considering D s psO.1 m 2 /s from our experimen- 
tal results, the average time the spins take to diffuse through the central suspended region is about 
x d = jjy ~ 30 P s f° r L=2.5 /im. Therefore we should see a decrease in the effective spin relaxation 
time Xfl t in case the average time the spins take to diffuse through the central region is in the order 
of the spin relaxation time of this region (t ; ), since most of the spins will relax before reaching the 
detector electrode. This consequence can be seen in the right graph of Figure 3c. If X{ is longer 
than x c i, the main limiting factor in the effective relaxation time is x . 

Now we will compare our theoretical results with the experimentally obtained values. As said 
before, we have a good representation of D s for the central suspended region from the values 
extracted by the experimental Hanle curves. The obtained values for the spin relaxation time, on 
the other hand, are apparently limited by the "dirty" outer regions. By performing Hanle precession 
measurements in fully LOR supported spin- valves, we obtain values of x s ps 150 ps and D s ps 0.02 
m 2 /s, which support our conclusions. This means that the spin relaxation time extracted from 
our Hanle precession measurements represents a lower bound on the spin relaxation time of the 
high quality central region. As far as we are aware this effect of inhomogeneity in a spin-valve 
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was never reported and it is of importance to experiments in any inhomogeneous system and not 
exclusively for graphene. For systematic studies in non-local spin valves as a function of the local 
properties of the central region (e.g. mobility, conductivity or carrier density) one must consider 
this inhomogeneity effects in the analysis, otherwise it may lead to erroneous interpretation of the 
experimental data. 

Although we are not able to determine the actual spin relaxation time in the high mobility 
region, the experimental conditions allow us to take a very important conclusion with regard to 
the spin relaxation mechanism. If we consider the D'Yakonov-Perel as the main mechanism for 
spin relaxation in graphene, we can assume that T s relates with the momentum relaxation time % p 
according to^ Equation 1. Since the charge diffusion coefficient relates to x p by: D c = T p , we 
can extract for our high mobility samples x p f^O.2 ps. We can then obtain an upper bound for the 
average spin-orbit coupling Aso <50 /leV using Equation 1. It is important to notice that the value 
for Aso that we estimate is an upper bound which is calculated from the experimentally determined 
lower bound on the spin relaxation time in the suspended high quality graphene. Applying the 
same procedure for typical SiC>2 supported graphene spin valves: T p =0.04 ps and Ts—200 ps we 
obtain an average SO coupling of Aso= 1 10 /leV. If this value was intrinsic for graphene we would 
obtain a strong reduction in the spin relaxation time in our high quality graphene devices down to 
Ti=50 ps. This value for t s is comparable to the time the spins take to diffuse through the central 
suspended region (t^) and smaller than the spin relaxation time for the non-suspended regions. 
This means that if the value for SO coupling obtained in the Si02 devices was an intrinsic value 
for graphene we should be able to observe a very low T s , which is not the case. From this result 
we can conclude that the nature of the SO fields in the graphene devices observed so-far is not due 
to graphene's properties, but due to extrinsic effects. The value obtained for Aso m ° ur suspended 
graphene devices is about a factor of two lower than what is usually obtained in the low mobility 
Si02 supported devices, although it is still five times higher than the theoretical limit of Aso ~10 
/ieV .I2J25J26] Q ur resu it s aj-e m agreement with recent theoretical predictions that Aso scales with 
the presence of adatomsJ-^S so higher spin lifetimes in high quality graphene flakes should be 
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expected. 

In conclusion we performed spin transport in suspended high mobility monolayer graphene 
devices contacted by ferromagnetic leads. We showed that electronic mobilities above 100,000 
cm 2 /Vs for the suspended region can be achieved in our devices via current annealing. Our mea- 
surements showed an increase up to one order of magnitude in the spin diffusion coefficient when 
compared to SiC>2 supported devices and very large spin relaxation lengths (A = 4.7 /im) were 
obtained. We did not observe a change in the measured spin relaxation time when comparing 
our results to the traditional Si02 supported devices. This effect is explained by considering a 
simple model that takes the inhomogeneity in our devices into account, considering not only the 
suspended regions of the devices, but also the supported part. We observe that the spin transport 
measurements of our high mobility graphene device strongly depends on how spins interact in the 
lower mobility graphene parts directly connected to it. In a similar way this inhomogeneity effect 
is expected to strongly affect the performance of a typical graphene nanodevice when connected to 
other graphene nanodevices or interconnects with different spintronic properties (due to possible 
edge roughness or other kind of spin scattering in the system). For future works, both in graphene 
and non-graphene based devices, one must take this effect into account in order to make concise 
conclusions. We were also able to give a higher bound for the spin-orbit coupling in our devices 
of 50 /leV, a factor of 2 lower than of those encountered in SiC>2 supported devices. 
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Supporting Information Available 
Device Preparation 

The devices are prepared in a similar fashion as the ones described by Tombros et al.,^ although a 
few changes were made in order to avoid degradation of our high resistive contact barriers. First a 1 
/im thick lift-off resist (LOR) film is spin-coated on a Si/Si02 (500 nm) substrate and the graphene 
flakes are exfoliated on top. Single layer flakes are then selected by optical contrast using a green 
filter.^122] p or ^ e i ec t ron beam lithography (EBL) process, to improve the undercut, we use a 
double layer Polymethyl methacrylate (PMMA) 50/410K resists dissolved in chlorobenzene and 
o-Xylene respectively. These solvents are used to prevent the removal of the LOR film during 
spin-coating. The contacts are then patterned and developed in n-Xylene (20° C). 

Using an electron beam evaporator with a base pressure lower than 8 x 10~ 7 Torr, we deposit 
0.4 nm of Aluminium followed by in-situ oxidation by pure Oxygen gas at a pressure higher than 
1 x 10 2 Torr for 15 minutes and the chamber is pumped down to the initial base pressure. This 
process is performed twice in order to get contact resistances higher than 10 k£L. After the high 
resistance barriers are deposited the chamber is pumped down to the initial base pressure and 60 
nm of Cobalt is evaporated. For some of the studied samples in this work the electrodes were 
capped by 3 nm of AI2O3 to prevent Co oxidation. The lift-off is done in hot (75° C) n-Xylene. 

To suspend the graphene flakes a second EBL step is performed with an area dose of 510 
/iC/cm 2 and developed in 1 -methyl 2-propanol. It was found that if the sample is immersed in 
Ethyl-lactat as described by Tombros et al.^ the AlOx barriers degrade, causing a very large 
increase in the contact resistance and loss of the spin-signal. After this final process the sample 
is bonded and loaded in a cryostat which is pumped down to a base pressure lower than 1 x 10 
Torr. 
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Current annealing 

After the sample is loaded in the cryostat and characterized, we perform a current annealing step 
to remove the impurities in the graphene flake and obtain a high mobility device. The whole 
procedure is carried at 4.2 K. 

To avoid degradation of the electrodes used for spin injection/detection, we apply the large 
DC bias for the current annealing in the two outer electrodes as depicted in Figure 4. The contact 
resistance of the inner contacts were measured before and after the current annealing step and 
showed no noticeable change. 




Figure 4: Top: A SEM picture of a typical device showing the schematics for the current annealing 
setup. On the left of the picture it can be seen two regions of suspended graphene flake, the 
outermost left successfully cleaned. On the right it can be seen two broken regions due to failing 
in the current annealing procedure. The scale bar measures 1 /im. Bottom: A cartoon illustrating 
the current annealing setup for comparison. 

We use a DC current bias-voltage compliance procedure to limit the power in our devices and 
avoid them to burn. The current is ramped up slowly (approx 1 jUA/s) until a determined value 
and then rapidly ramped down, at a rate 4 times faster than the ramping up. After each sweep in 
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current we check the gate-voltage dependence on the sheet resistance to keep track of the device's 
mobility. For our samples, this current annealing procedure had a success rate of about 33% (4 out 
of 12 regions showed high mobility), comparable with our previous results.^ One of the regions 
was current annealed twice. The first procedure resulted in a mobility of /i ~ 10 5 cm 2 /Vs, and 
after the second current annealing the mobility was improved to /! « 3 x 10 5 cm 2 /Vs. Despite 
performing spin-transport measurements in this sample, we also obtained spin signals in another 
sample with similar properties. 



Details on the simulation 

To represent our sample we extended the ID model by Popinciuc et al. 7 to include three different 
regions: two semi- infinite outer parts of width W sandwiching one inner part of length L and 
width W. In this model we solve the stationary Bloch equations for the three components of the 
spin accumulation ]l s (x) in the presence of a magnetic field B: 

D s V%-^ + rBxfl s = (2) 

where D s is the spin diffusion constant, is the spin relaxation time and y = gji^t i~ [ is the gyro- 
magnetic ratio. All the parameters of the equation above can be set separately for each of the three 
regions, but for simplicity we make the two outer regions identical. For the boundary conditions 
we take: 

fi s (x = ±oo) =0 

fi s (x = + ) =/i,(x = 0_) 

H s (x = L+) = n s (x = L_) 

O _ G W d^l s (x=0-) _ (jjW dfl s (x=0 + ) . /Xs(x=0) 
" 2e dx 2e dx 2eR c \ 

n _ QjW d^l s (x=L-) _ q W d^i s (x=L + ) fi s {x=L) 
2e dx 2e dx 2eR c 2 
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where j8 = ^ for the x component of ji s and /3 = for the y and z components, with P being the 
spin polarization of the charge current / injected in the left boundary. The contact resistance of 
the contact at the left (right) boundary is represented by R c \( c 2), and the conductivity of the inner 
(outer) regions by <% ). The contact induced spin relaxation is represented by the last term of items 
4 and 5,^ although contact effects were found to be negligible for our results when we consider 
values obtained experimentally in our samples. 

We have to study the effect of four different parameters: the spin diffusion constants A and 
D a , and the spin relaxation times T, and % , where the subscripts "i" and "o" refer to the inner and 
outer regions respectively. In order to be able to observe the effect of each one of the parameters 
separately we calculated several Hanle precession curves keeping three of them constant and vary 
the remaining one. The simulated precession curves for a few sets of parameters are depicted in 
Figure 5. These curves were then fitted using the solution for the Bloch equations in a homoge- 
neous system, like we fit our experimental results. From these fits we obtain an "effective" spin 
diffusion constant Dfj t and relaxation time as shown in the main manuscript. We also tried to 
fit our experimental data with the curves we get from our model, but it lead to similar results to 
the ones obtained using the solution for a homogeneous system. It is worth noting that when we 
compare the curves in Figure 5a, we observe that in the case of changing D the obtained preces- 
sion curves change in magnitude but not in shape. This means that we do not observe any change 
in the values obtained for Dfj t or Ty#. On the other hand, analyzing Figure 5b we see that changes 
in Dj does not only change the magnitude of the spin signal but also the shape of the curve, which 
results in changes in Dfn with changes in £>,. A similar effect is observed in the analysis of Figure 
5c and Figure 5d. 



The effect of the sheet resistance on the Hanle precession 

The same way as we can change the values for the spin diffusion coefficients and relaxation times 
for the inner and outer parts (D a , Du T and T,) we can also change the values for the square 
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Figure 5: The calculated Hanle precession curves for different values of (a)D , (b)D,-, (c)T and 
(d)r,, while keeping the other parameters fixed. 
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resistances of the inner and outer regions, Rj and R a respectively. By applying the same procedure 
of generating a Hanle precession curve and fitting it with the homogeneous model we can extract 
the effective spin relaxation time xy, ? and the effective spin diffusion coefficient Dfu- The results 
for different combinations of R\ and R G are presented in Figure 6 below. 




C/3 

N 

E 



C/5 
CL 




Figure 6: (a) The spin signal obtained for different values of the inner (outer) sheet resistance while 
keeping the outer (inner) constant at 5 \SL. The results obtained for Dfn and xy ;f for changing (b)i? ; 
and (c)R a . The outer or inner sheet resistance was kept at 5 kQ. while the other was changed. The 
values for the diffusion coefficients and spin relaxation times are: D o =0.01 m 2 /s, D,-0.1 m 2 /s, x () = 
200 ps and T,= 1 ps. 



As it can be seen in Figure 6 (b) and (c), the conclusions obtained in our main text that Dfn 
is determined mainly by the inner region and Ty# by the outer region remain unchanged when we 
consider resistances in the range of those we encounter experimentally (1 to 5 k£l). 
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